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ABSTRACT 


This report describes a newly developed non-isotropic multiple-scale turbulence model 
(MS/ASM) for complex flow calculations. This model focuses on the direct modeling of 
Reynolds Stresses and utilizes split-spectrum concepts for modeling multiple-scale effects 
in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows 
show that the current model perform significantly better than the single-scale k - e model. 
The present model is relatively inexpensive in terms of CPU time which makes it suitable 
for broad engineering flow applications. 

INTRODUCTION 

Presently the great majority of engineering flow computations are based on an eddy 
viscosity linkage through Reynolds-stresses/Mean-strain rate constitutive relationship, the 
most often used version being one in which the eddy viscosity is obtained from local 
values of the turbulence kinetic energy (k) and its dissipation rate (e) [1], The k and e 
are obtained from transport differential equations solved simultaneously with the mean 
momentum equations. Some of the shortcomings inherent in this type of k - e models are 
: (a) the use of one set of characteristic scales to represent a spectrum of eddy sizes and, 
(b) the isotropic assumption of the gradient type eddy viscosity Reynolds-stress/strain 
rate constitutive formulation. Experimental observations indicate that energy-containing, 
turbulence generating eddies are larger and exhibit a different rate of development than 
do smaller turbulence energy dissipating eddies. In the classical turbulence models such 
as the k - c model and the Reynolds stress models [2], only the generation and dissipation 
of turbulent kinetic energy are considered while the cascading processes are assumed to be 
in equilibrium. 

Besides, the effect of rotation on the turbulence, especially through the modifica- 
tion of the turbulent energy cascading process, has not been explored by the conventional 
single-scale models. Experimental data shows that the predominant effect of rotation is 
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to decrease the rate of dissipation of the turbulence and increase the lengthscales, espe- 
cially those along the axis of rotation. Due to the generation of inertial waves, the net 
energy transfer from large eddies to small ones was reduced. To account for the different 
mechanisms involved in the evolutions of different eddies, the generation, cascade, and 
dissipation of turbulent kinetic energy have to be considered. 

In this study, we have developed a multiple-scale turbulence model based on the 
split-spectrum methodology of Hanjalic et al [3]. This model is derived by partitioning 
the turbulent kinetic energy spectrum into a set of wave number regions, by integrating 
the partitioned energy spectral density equation to obtain the energy evolution for each 
wave number region, and by considering the energy transfer rate between the adjacent 
wave number regions. Recent developments including the rotation effect and anisotropy 
effect will be described. Model validations including homogenous rotating flows, free shear 
flows and recirculating flows will be compared with available experimental data and other 
single-scale turbulence models. 


THE TURBULENCE MODEL 


The underlying idea of the multiple-scale model has been described in [3,4]. To sum- 
marize, the closures and the convection-diffusion equations for the model are described 
below. The turbulent kinetic energy and the energy transfer rate equations for the energy- 
containing large eddies are given as 
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dt 

dt 
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where P r is the energy production rate represented by P r = — TTiujdUi/dxj . The 
turbulent kinetic energy equation and the dissipation rate equation for the fine scale eddies 


are given as 
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In the above equation, fl is the local mean vorticity or angular velocity of the rotating 
frame. The model constants in equations (1) - (4) have been determined by transform- 
ing the modeled turbulence equations into asymptotic turbulence decay rate equations 
of grid-generated homogeneous turbulence and other simplified flow situations. Detailed 
derivation of the system of equations for the model constants which is similar to that in 
Ref. [4] can be found in Appendix A. These model constants are given in Table 1. 

If the above turbulence model equations are used to construct an eddy viscosity based 
on Boussinesq assumption, the eddy viscosity formula proposed here is given as 

= c- ( 6 ) 


and the Reynolds stress tensor is expressed as 


v~uj = Rij — -vtSij + 


(7) 


in which 


Sii - 


l,dUi , dUj 


+ 


2 dxj dxi 


k = -uiui 
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( 8 ) 

(9) 


As discussed widely in the literature, the linear isotropic constitutive relationship 
used in equation (7) yields highly inaccurate predictions for the normal Reynolds stress 
differences. For example, in a fully developed turbulent channel flow, equation (7) predicts 
that the normal Reynolds stresses are all equal, i.e., that 


Rn = R22 = R33 («i«i = u 2 U2 - 113U3) 


which indicates isotropic turbulence and is in contradiction to experimental observation. 
In this study, the non-isotropic aspect of the turbulence will be modeled based on simpli- 
fication of the second-moment closure methodologies. In the second- moment closure, the 
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transport equation for Reynolds stress components R,j was derived [5,6] : 

§- t {pRa) - By = Pa + ~ “i (>°) 

in which the terms on the left hand side of (10) represent the convective and diffusive trans- 
port of Rij, respectively, and the terms on the right hand side represent the production, 
pressure-strain (redistribution), and viscous dissipation tensors of Reynolds stress. Among 
these terms, closure assumptions are required for Dij, $ij and e^. To date, equation (10) 
still represents the most complex turbulence models used in engineering calculations for 
which the non-isotropy and extra complex strains can be automatically accounted for. 
However, the extra Reynolds stress differential transport partial differential equations (six 
of them for 3-D flows) required to be solved introduce severe computer cost penalty. As 
a result, the approach used in this study to retain the advantages of this level of model- 
ing is to simplify the differential transport equations such that they reduce to algebraic 
expressions, the so called Algebraic Stress Model (ASM)[7j. 


There are several approaches for algebraic stress model formulations. The one adopted 
here is the similitude analysis proposed by Mellor and Yamada [ 8]. Their formulation leads 
to linear algebraic equations for the stresses based on a rigorous similitude analysis of the 
transport equation for anisotropy (uiu] - 2/36ijk). This approach does not subject to the 
limitation that uiuj/k = constant assumption invoked by Rodi [2]. Equation (10) was 
thus simplified to : 

P i} - hijPk + $ij 3S 0 (11) 


It should be noted that the model for should be consistent with the order of magnitude 
analysis which leads to equation (11). The model proposed in this study has the following 
form : 

$ij = -Cip^-(u7u7- ^Sijikp + k t )) - c 2 {P t j - -t>ijPk) (12) 

and 


Pij = Ui Uj 


dUj 

dxj 


(13) 
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(14) 


P * = \ P " 

The first term on the right hand side of equation (12) contains only turbulent velocities 
and is formulated after Rotta’s linear return-to-isotropy model [6]. Note that the time scale 
used in the turbulence time scale associated with the high frequency eddies in the multiple- 
scale model. Detailed derivations can be found in Appendix B. The above model causes 
a tendency towards isotropy if c\ is greater than unity. In this study, Ci is chosen to be 
1.8 - 0.5 / in which / is a wall correction term affecting the normal redistribution actions 
of the pressure-strain correlation. This ’’echo” effect has been discussed extensively in the 
literature [5,6]. We use the following simple proximation function : 

/ = !-(£)’ 

where y is the distance from wall and D is the half width of the duct passage. 

The second term on the right hand side of equation (12) is referred to as the ’’mean- 
strain” part of the pressure-strain tensor. We model this term based on the ’’isotropization 
of production” of Ref. [5]. As seen from the formulation, this term associates with the 
production of turbulent kinetic energy and thus links with the time scale of the energy 
containing eddies. The constant used here is c-i = 0.85 — 0.06/. 

For wall bounded flows, a ratio of k t /k p was derived based on the wall function ap- 
proach [2]. For practical purposes, these wall functions were derived from the near- wall 
equilibrium assumptions and logarithmic velocity profiles. Detailed derivations and for- 
mulations can be found in Appendix C. 

RESULTS AND DISCUSSIONS 

For model validations, several testing cases including grid-generated turbulence with 
rotation, free shear flows and two-dimensional as well as axi-symmetric recirculating flows 
are calculated. 
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Grid- Generated Rotating Flows 


The structure of the turbulence is changed through the action of the Coriolis forces : 
the turbulent length scales are increased in the rotation axis and the energy cascading from 
larger toward smaller dissipative eddies is inhibited [9]. Figure 1 compares the numeri- 
cal predictions of the multiple-scale model with corresponding experimental conditions of 
Wigland and Nagib [10] as listed in Table 2. The model provides quite accurate predictions 
of the decay of turbulence intensity at different rotation rates and different flow conditions. 
The agreement is good for all cases, even when in case (2) the turbulence decays faster for 
9, = 20s -1 due to different initial time scales. Under the same initial conditions, the model 
predicts slower decay of the turbulence with increasing rates of rotation, in agreement with 
the data. It should be noted that the standard k — e model is insensitive to changes in 
rotation rates because it has no provision to account for this effect. 

Free Shear Flows - Jets 

The necessary step in turbulence model development is to use the model to predict 
free shear flows. Here, three different jet flows in the fully developed region (self-similar 
regions) are calculated using an accurate finite different scheme based on Keller’s Box 
method. The set of governing partial differential equations were transformed to a set of 
coupled ordinary differential equations. Such as the final coupled ordinary differential 
equations for the round jet is shown in Appendix D. Details of the numerical scheme and 
the method of calculating sharp boundaries between turbulent shear flows and irrotational 
ambient fluid can be found in [11]. The multiple-scale model used here is the one associated 
with an eddy viscosity formulation of equation (6) (the MS/EV model). 

The computed values of the spreading rates of the plane jet, round jet and radial jet 
using the MS/EV model are given in Table 3. Also shown for comparison are the computed 
spreading rate using various single-scale models as well as the experimental values for these 
quantities. The spreading rate of the round jet is experimentally observed to be about 20 
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% less than that of the plane jet. The inclusion of the rotation time scale in the MS/EV 
model has the effect of increasing the energy transfer rate from energy-containing eddies 
down to the dissipative eddies due to the vortex stretching mechanism and produced less 
spreading for the round jet. The MS/EV model gives the best overall predictions in terms 
of spreading rate for these jet flows. 


In Figure 2-4, the dimensionless mean velocity of the three jets are plotted versus the 
similarity coordinates in which the transverse locations are normalized with the velocity 
half-width. The mean velocities calculated using various models are all in reasonable 
agreement with the data over the central position of the jet but deviates towards the edge of 
the jets. This is due to the intermittence close to the free turbulent/non-turbulent interface 
and the turbulence models do not model these unsteadiness explicitly. The MS/EV model 
appears to give better predictions than the others and gives a little wider jets. 

2D Backward-Facing Step Flow 

For this and the following case, the governing partial differential equations were solved 
in the fully elliptic form. The differential equations were discretized based on the control- 
volume approach with a non-staggered grid arrangement. A second order upwind dif- 
ferencing scheme is employed for convection terms. A time-marching sequence employs 
the PISO[12] algorithm to resolve the momentum equation using the pressure incremental 
equation which enforces continuity at each time step. Calculations are carried out until 
steady state solutions are reached. Details of the numerical method can be found in [13]. 

The predicted values of the reattachment length of the recirculation zone behind the 
step were made using the standard k - e model, the single- scale Algebraic Stress Model 
{k-e / ASM), the multiple-scale model with eddy viscosity formulation (MS/EV) and the 
most recent non-isotropic multiple-scale model (MS/ASM) are listed in Table 4. It can be 
seen that the single scale models (k-e and k-e /ASM) underpredict the reattachment 
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length while both the multi-scale models give good predictions of the reattachment length 
compared to the experimental data. 

In terms of turbulence quantities, comparisons of the predicted streamwise and lateral 
turbulence intensities as well as the Reynolds stress component uv obtained from both the 
k — e /ASM and the MS/ASM are shown in Figure 5. The results were normalized with the 
inlet centerline velocity. Non-isotropic turbulence structures are predicted by both models 
with the MS/ ASM doing a better job with exceptions of near- wall regions. The inability 
of the models to reproduce the near-wall turbulence correctly behind the step is traced to 
the weakness of the wall proximate term of the pressure-strain model. 

Axisvmmetric Pipe Expansion Flow 

The 1:3 pipe expansion flow of [14] is selected for model validations. Again, the 
reattachinent length predictions from the four models are listed in Table 5. Both the 
MS models give satisfactory results. The MS/ASM actually predicted a slightly shorter 
reattachment length in comparison with the MS/EV model. The results of the k - e/ASM 
and the MS/ASM predictions of streamwise and radial turbulence intensities are illustrated 
in Figure 6. It is interesting to note that the strength of longitudinal and radial intensities 
(just after the pipe expansion x/H = 2 ) are very well predicted by the MS/ASM model. 
Within the recirculation regions, both models gave reasonable features of nonisotropic 
turbulence with the MS/ASM doing better performance. After the flow reattachment, the 
models have the tendency to predict the ” return to isotropy” faster than the experiments. 
In the MS/ASM model, although a different turbulence time-scale is used to model this 
mechanism, the coefficient is not returned, i.e., the coefficient used in the conventional 
Rotta model is adopted without change. This aspect of modeling requires further research. 
The predicted turbulent shear stress from the two non-isotropic models are compared with 
the measurements in Figure 7. Overall, the MS/ ASM model does a better job of predicting 
the Reynolds stress component. 
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CONCLUSIONS 


Validation studies of the newly developed turbulence model have been presented for 
rotating flows, free shear flows and recirculating flows. It is evident that in general the 
MS models perform better than various single-scale models. The MS/ASM model also 
shows good capability of giving quite accurate predictions even in second order turbulence 
quantities for the engineering calculations. Besides, the CPU time and Computer storage 
for the present model are relatively inexpensive (about 20 % more for 2-D backward- 
facing step flows) compared to the k — e model. Some discrepancies between the data 
and the MS/ASM solutions can be related to the pressure-strain models, especially the 
wall-damping model. Further validations are required to further develop the model. 
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NOMENCLATURE 

D jet exit diameter or width 

Dij diffusion term of Reynolds stress flow 

/ wall proximation function 

H step height in backward-facing step flow or expansion 
k p turbulence kinetic energy in the production range 

kt turbulence kinetic energy in the dissipation range 

Pij production term of Reynolds stress 

P r kinetic energy production rate 

Rij Reynolds stress tensor 

U c mean centerline velocity of jets 

U j Xj component of mean velocity 
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Uj 


y 

h 


e p 


Ct 


1 


Vt 


VK , &k t 






Xj component of velocity fluctuation 

coordinate direction, 1 - streamwise, 2 - transverse 

streamwise coordinate direction 

transverse coordinate direction 

kronecker delta function 

viscous dissipation tendor 

transfer rate of kinetic energy 

dissipation rate of kinetic energy 

similarity coordinates 

turbulent viscosity 

constants in model equations 

constants in model equations 

pressure-strain tensor 
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APPENDIX A. THE TURBULENCE MODEL CONSTANTS 


Determination of the turbulence model constants is based on the assumption that the 
ratio kt/k p will asymptotically approach constant values in simple homogeneous turbulent 
flows such as the decay of grid turbulence, a nearly homogeneous shear flow, and the near 
wall equilibrium turbulent flows. The model constants also need to satisfy the realizability 
conditions (i.e., the turbulent kinetic energies, energy transfer rate, and dissipation rate 
cannot become negative). This procedure is described as following. 


At hight Reynolds number, the multiple- scale model equations (l)-(4) for the homo- 
geneous flow will be deduced as follows: 


— (7n,fle r 
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The above set of equations can be solved analytically. The final results 
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For the grid turbulent flow without rotation, the results become 

r L = < 1 + f r ' 

KpQ tQ 


r = ( 1 + 7) 

mo r o 

it- = (i + !)-(«+» 

CpO <0 
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where 


t n — 7T 


0 

e p0 


w ^to 

t 0 = m — 

f to 


/ 0 and /q stand for the time scale of large eddies and small eddies, respectively. It can 
be seen from equation (A-9) that the value of n is the decay rate of kinetic energy k p . 
According to the experimental data, n = 1.11 ~ 1.25, hence C P2 = 1.8 ~ 1.9. To account 
for the characteristics of the energy spectrum kt/k p , C P2 is set to be the function 1.90 *(1 — 
0,2 ■ ktjk p + 0.2 * fcjV&p), In this function, the value of C P2 is always within the experimental 
range if k t /k p is less or equal to one. The functional form of C P2 is equivalent to t he variable 
energy transfer function proposed by Hanjelic et al. 


To recover the k — e model, the constant C pi is given as conventional value 1.42. If the 
partition is moved to too high wave numbers, i.e., kt = 0, then the multiple-scale model 
may reduce to k — e model. 


In order to keep the same decay rate of e p as that of e*, we set C P3 = C t3 . The value 
of C P3 and Ct 3 is determined as 0.042 by the computational results to fit the experimental 
data for homogeneous rotating flow. In this study, Ct 3 is chosen as a function of 0.042 • 
(1 + 0.08 -€t/€ p ). The ratio of e* and e p is a characteristics for degree of spectral imbalance. 
Inclusion of t t /e p into the functional form of Ct 2 (and Ct x ) allows the ’’inter-talk” between 
the large eddies and small eddies. 

The values of C tl and C tl were optimized to give good agreement of the computational 
results with the experimental data for the free shear flows and 2D recirculating flows. 
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According to the observation the effect of rotation on the turbulence is to decrease the 
rate of dissipation of the turbulence and increase the lengthscale, especially those along 
the axis of rotation. Due to the generation of inertial waves, the energy transfer from large 
eddies to small eddies is reduced. Then the additional time scale to account for the effect of 
rotation is introduced in the equations of e t and e p , respectively, fi is the angular velocity 
of the rotating frame and the local mean vorticity, that is defined as following form: 

n = (n«n<i) I/2 , = 

Finally, it is mentioned that the prandtl numbers for turbulence quantities are set 
following [4], they are 

Vk,, = = 0 . 75 , < t €i , = <r ft = 1.00 
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APPENDIX B. THE ALGEBRAIC STRESS MODEL FORMULATIONS 

The transport equation for Reynolds stress components Rij is given as 

£- t (pRij) - o tj = Pa + - * ,j (b- i) 

in which the terms on the left hand side of equation(B-l) represent the convective and 
diffusive transport of Rij, respectively, and the terms on the right hand side represent the 
production, pressure-strain and viscous dissipation tensors of Reynolds stress. Using the 
similitude analysis proposed by Mellor and Yamada [8], equation (B-l) can be simplified 
to 

Pij-pijPk + ^ij=0 ( B-2 ) 

The model presented in this study has the following form: 

j = — ciPT-( u i u j ~ ~6ij(k p -(- k t )) — c 2 {Pij — ~f> ijPk ) (B ~ 3) 

kt O ‘J 

and 

(B-4) 

Pk=\Pu (B-V 

Then equation (B-2) can be written as 

Pij - hijPk - c lP € ±{u^-- \Sij{k p + *«)) - c 2 (P tj - hijPk) = 0 (B - 6) 

or 

(1 - c 2 )Pij -I- \Sij(c 2 Pk + c lP j-{k v -I- kt) - Pk) = c ip ^-uiu] [B - 7) 

i Kt Kt 

For 2D turbulent flow, final Reynolds stress components u\, it?, and il x u 2 were ob- 
tained as following: 

— — 2k[(p/3 • (ci • e</ €p 4- c 2 • kt/k kt/k) fcf/fe • (1 C 2 )^i lt 2 ‘ dU\/dx 2 \ , p 

“i = " c 1 e 1 + 2Ml-C2)-5U 1 /ax 1 
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— _ 2fc[cp/3 ■ (ci • e t /e p + cz • k t /k — k t /k) — k t /k • (1 — c 2 )uiu 2 ■ dU 2 /dx\] 
2 c i € t + 2k t (l — c 2 ) • dU 2 j dx 2 


(B- 9) 


UiU 2 = 


Cl ft 


-dV 2 -dlh. 

( C 2 - l){u(— + u- 2 — ) 


dx 2 


(5-10) 


where 


k — kp 4 kf 
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APPENDIX C. WALL FUNCTION BOUNDARY CONDITIONS 


To avoid the need for detailed calculations in the nearwall regions, equations are 
introduced to link velocities, & p , k t , e p , and on the wall to those in the logarithmic region. 
In this region, the flow is dominated by a one-dimensional Couette flow characterization of 
the flow diffusion perpendicular to the wall. Also it provides a way around this region of 
steep nonlinear and variation of the variable and the fact that laminar and turbulent effects 
become of the same order of magnitude. The new equations introduced in the momentum 
equations and the turbulence model equations are used in the finite difference calculations 
at nearwall points. 

(1). The universal velocity profile 


K = b n <^ + ) (C- 1 * 

where 

Vr = C-j)' 2 (C- 2 ) 

U r is the friction velocity; r w is the wall shear stress; = yu r ju is the non-dimensional 
wall coordinate; k is the Von Karman constant, k — 0.4187; E is a experimentally deter- 
mined constant coefficient, E — 9.0. 

(2). The wall shearing stress 

Assuming r* is an approximation for r w very near the wall. The r*. is formulated by 
observating that convection and diffusion of turbulence kinetic energy are nearly always 
negligible in this region. Deleting these terms from the ^-transport equation and invoking 
isotropic viscosity leads to 

n = Cl /7 pk (C — 3) 

From equations (C-l) and (C-2), thus we obtained 

= -[ P KC 1 /*k l ' 2 /]n{Ey + )] -U (C- 4) 
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where the negative sign is inserted since t w and U must have opposite directions. 

(3). The total kinetic energy 

From equation (C-2), the total kinetic energy very near the wall is given as 

k = C~ l/2 Ul (C — 5) 


(4). The energy dissipation 


As for energy dissipation near a wail, the length scale is assumed to be proportional 
to the normal distance from the wall, that is 

A; 3 / 2 


= Ci-y 


(C- 6) 


which leads to 

€ = — C 3 J 4 k 3 / 2 ( c - 7) 

Ky ^ 

This is the effective wall boundary condition on e. 

(5). The kinetic energies k p and k t 


k = k p - (- k t 


(C- 8 ) 


If deleting convection and diffusion of turbulence kinetic energy and assuming P r = 
e p = e, in the region very near the wall , from equations (2) and (4) the ratio of k t 
and k p is derived 


kt 

Vp 


C tl -C tl -C t3 Clk t /e t 

Cpi — C P7 — Op 3 Clk p / e p 


(C- 9) 


when fl ~ 0, 


kt C tl - C t 


k r 


- C' n 


( c - 10 ) 


tp V- p! 2 

Equations (C-8) and (C-9) or (C-10) provide two constraint conditions for the kinetic 
energies k v and k t . 
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(6). The Transfer rate e p and dissipation rate e t 

In the near wall region, where the turbulence is in equilibrium (the production rate is 
approximately equal to the dissipation rate), it can thus be assumed that 

Pr = €p = e t (C — 11) 
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APPENDIX D. SIMILARITY EQUATIONS FOR THE ROUND JET 


The free jet flow becomes self-preserving at some distance away from its origin, sim- 
ilarity solutions for mean quantities and turbulence quantities exist. By proper choice of 
local length scale and velocity scale, the set of momentum and turbulence quantity equa- 
tions can be transformed to a set of coupled non-linear ordinary differential equations. 
The resulting equations are quite lengthy and only the case of a round jet is summarized 
below. 


Similarity equations for the round jet 


\ flf , r/f \'l 1 I 

— if — + /(— ) - ~ T u =r u 
77 77 77 77 

-“(2 f'p + fp') =(—)'r u - m + -T fcp + T k 
77 V V 

~ “(2 f'q + fq) =m - n + -r kt + r kt 
77 77 


kp 


(D- 1 ) 
(D- 2) 
(D- 3 ) 


--( 4 /'m + fm') =C p it* — (— )' - C p2 { 1 - 0.2- + 0.2^-) — 

7] p X T] / p P P 

- C p 3 m[X 2 V f" + (-)'] + -r ep (D - 4 ) 


77 77 

2 


--( 4 f'n + fn') =(7(1 ( 1 -1- 0.01 — )— - C i2 — 

rf m q q 


- C t3 ( 1 + 0 . 08 -)n[A 2 77 /" + (£)'] 4- — T et (D - 5 ) 
m r/ 1 ] 

A' =0 (D-6) 


Boundary conditions: 


/' 


at 77 = 0 : — = 1; / = p = q' = m' — n = f" = 0 

V 

at 77 = 1: TTl =0; P ~ q = TTl = 7Z = 0 


(D-7) 

(D- 8) 


Definition of variables: 


()' = /(®i) = A®i, i 7( = C m ^-^- 

dr] m 
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k p _ kt_ 

p- uy q ~ ir- ’ 


m = 


u 3 ’ 


n = 


e </(x i) 

U} ’ 


77 = 

iTo 

lh 

t 

Tu = V t { — )\ 

n 


Uc 

n ’ 

Vt , 
Tkp = P > 

Ckp 

Tkt = 

Vt , 
9 1 

&kt 

T *P 

Vt , Vt 

= m , T (t — — n 

<T € p 


It would suffice to say that due to the complete neglect of molecular diffusion, the 
choice of the characteristic length scale corresponds to the locus of points where the tur- 
bulent diffusivities are zero. The location of these points must be calculated as part of the 
solution to the problem. The existence of this sharp boundary is insured by requiring the 
turbulent diffusivity {v t ) to be zero there. 
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TABLE 1 


Constants and Functions of Multiple-Scale Model 


Constant or Function 

Value or form 

Cpt 

1.42 

Cp, 

1.90(1 — 0.2ki/k p + 0.2fc( /fcp) 

CP3 

0.042 

Ct x 

0.96(1 +0.01t t /e p ) 

Cu 

1.12 

Ct, 

0.042(1 +■ 0.08c t /e p ) 

*k v 

0.75 

<?k, 

0.75 


1.00 


1.00 




TABLE 2 


Experimental data of Wigeland and Nagib (1978) at initial 
time (Ut/m=20) and parameter ranges covered 




Ca.se 

(1): 

M 

= 0.0039 (m) 


n 

(-S’ 1 ) 


0 


20 

80 

u 

(m/s) 


8.69 


8.66 

8.82 

u. 

(m/s) 


8.69 


8.67 

8.96 

k 

(m 2 /s 2 ) 


0.10625 


0.10105 

0.1157 

e 

( m-/s 3 ) 


17.67 


16.45 

18.19 



Case 

(2): 

M 

= 0.00625 (m) 


n 

(^I 1 ) 


0 


20 

80 

u 

(m/s) 


5.33 


5.45 

5.48 

U, 

(m/s) 


5.33 


5.46 

5.71 

k 

(m 2 /s 2 ) 


0.0425 


0.0453 

0.04975 

e 

( m 2 /s 3 ) 


2.649 


3.591 

3.300 



Ca.se 

(3): 

M 

= 0.00254 (m) 


n 

(s' 1 ) 


0 


20 

80 

u 

(m/s) 


8.51 


8.58 

8.79 

Us 

(m/s) 


8.51 


8.59 

8.93 

k 

(m 2 / 5 2 ) 


0.0885 


0.0887 

0.09845 

€ 

(m 2 /, 3 ) 


24.27 


22.56 

22.42 


TABLE 3 


Spreading Rate Predictions Versus Experiments 


Model 

Plane Jet 

Round Jet 

Radial Jet 

k - e [3] 

0.110 

0.125 

0.094 

Morese [5] 

0.102 

0.085 

0.099 

Pope [2] 

0.110 

0.091 

0.042 

Hanjalic & Launder [6] 

0.119 

0.107 

0.077 

Chen & Guo 

0.111 

0.084 

0.094 


Experiment 


0.100 - 0.110 


0.086 - 0.09 


0.098 - 0.110 


TABLE 4 


Reattacliment length for Backward facing step turbulent flows 
9:1 Experiment Driver and Seegmiller(1986) 


Experiment 

k — e Multiple-Scale 

k - tj ASM 

MS / ASM 

- 6.1 

4.76 5.94 

4.65 

5.70 

TABLE 5 

Reattachment length for Axissymatrical Recirculating turbulent flows 


3:1 (Axis.) Experiment 

Najed(1987) 


Experiment 

k — e Multi pie- Scale 

k - c/ ASM 

MS / ASM 


~ 9 


7.10 


8.82 


6.80 


8.70 







model [ 15 ] 
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omparision of computational and experimental 


. Pope [17] 
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Figure. 3 Comparision of computational and experimental mean 
velocity for the round jet 



Figure. 4 Comparision of computational and experimental 
velocity for the radial jet 



Figure. 5 Reynolds stresses profiles for the backward facing step 
turbulent flow (9:1) 



Figure. 6 Axial and radial turbulence intensities profiles for the 
recirculating turbulent flow (3:1, Axis.) 
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